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Abstract Bodies of density one half (of the fluid in which they are immersed) 
that can float in all orientations are investigated. It is shown that expansions 
starting from and deforming the (hyper) sphere are possible in arbitrary dimen- 
sions and allow for a large manifold of solutions: One may either (i) expand 
r(n) + r(— n) in powers of a given difference r(u) — r{— u), (r(n) denoting the 
distance from the origin in direction n). Or (ii) the envelope of the water planes 
(for fixed body and varying direction of gravitation) may be given. Equivalently 
r(n) can be expanded in powers of the distance h(u) of the water planes from 
the origin perpendicular to u. 

1 Introduction and Results 

A long standing problem asked by Stanislaw Ulam in the Scottish BooklJ (prob- 
lem 19) is, whether a sphere is the only solid of uniform density which will float 
in water in any position. Such a solid is called a floating body of equilibrium. It 
will be in indifferent equilibrium in all orientations. 

The simpler, two-dimensional, problem to find non circular cross-sections of 
a long cylindrical log which floats without tending to rotate (the axis of the log 
is assumed to be parallel to the water surface.) was solved for relative density 
Pd = 1/2 in 1938 by Auerbach [2]. He found a large class of solutions. This is 
in contrast to the solutions found for pd 7^ 1/2 in two dimensions [U 01 [5] and 
for central symmetric bodies in three dimensions [6], where the variety of shapes 
is much more restricted, provided one restricts to star-shaped bodies. 

Here it will be shown that also in dimensions d > 2 there is a large variety 
of bodies which can float in arbitrary orientation at pd = 1/2. More precisely: 
Denote the distance from the origin to the boundary of the body by r(n), with 
the unit vector n pointing into the direction. Starting from the hypersphere 
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r(n) = ro a deformation governed by an expansion parameter e will be consid- 
ered 

oo 

r(n) = r + ^e fe r fe (n). (1) 
fc=i 

Due to Archimedes' principle one half of the body is below the water, the other 
half above. Thus the water plane cuts the body into two halves. Denote the 
height of this plane above the origin by h(n) with n perpendicular to the plane, 
then h(— n) = — h(n). A similar expansion will be performed for /i(n), 

oo 

h(n)=J2^hk(n). (2) 

fc=i 

For practical reasons restriction will be made to solutions which obey 

rfc(-n) = R fc r fc (n) (3) 
/i 2 fc+i(-n) = -h 2k+1 (n), h 2k {n) = 0. (4) 

In each order in e there are two conditions on and h/.. One condition guaran- 
tees that the volumes above and below the water level are equal (V-condition, 
cq. ITTj) . the other guarantees that the potential energy does not depend on the 
orientation (Z-condition, eq. I70p . 

With the restrictions (01 2]) one finds in odd orders e 2fc+1 that the Z-condition 
is identically fulfilled. The V-condition yields that a linear combination of r 2 k+i 
and h 2 k+i is determined by r n and h n with n < 2k + 1. In even orders e 2k the 
V-condition is identically fulfilled. The Z-condition determines r 2k . 

This may be used in the following ways: 

Given r. If ri(n) is given, which obeys n(— n) = — ri(n) and r 2 k+\ = 
for k > 0, then the conditions yield in odd order hk and in even order r^. If 
one chooses also r\ = 0, then the r 2 k will also vanish (except n independent 
contributions). This is in agreement with the theorem by Schneider pj [8] and 
Falconer [9], also referred to by Hensley in the Scottish book lj, compare [6j: For 
arbitrary dimension d and density 1/2, if the body is star-shaped, symmetric, 
bounded and measurable, then it differs from a ball by a set of measure 0. 

Given h. If fta(n) is given, which obeys hi(— n) = — fta(n) and hk = for 
k > 1, then the conditions yield ru in all orders. Suppose the envelope of the 
water planes is given. The distance of the planes tangent to this envelope yield 
ft,(n). Thus to a given envelope one may construct the surface of the floating 
body. If one rescales r and h by a factor 1/e, then one obtains 

oo 

r(n) = ^+^-V fe (n), h(n)=h 1 (n), (5) 

6 k=l 

Thus for a given envelope one obtains a one-parameter manifold parametrized 
by e. These envelopes have the property that there is exactly one tangent hy- 
perplane to them parallel to any given plane. As a consequence these envelopes 
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have cusps and are quite different from surfaces of convex bodies, which have 
two tangent hyperplanes parallel to any given plane. A one-parameter mani- 
fold for a given envelope was obtained in dimension d — 2 by Auerbach[5]; see 
appendix lA.il From given h(n) one can determine the envelope fappendix lA.2[) . 

Basically the expansion scheme of ref. [B] will be used. In subsection 12 . 1 1 the 
two conditions for equilibrium is given. Two expansions are used: the expansion 
in e, which describes the deviation from the hypersphere (]2.2[) and the expansion 
in ultraspherical harmonics (|2.3[) . In subsection l2.4l thc expansion of V and Z in 
ultrasphcrical harmonics is completed. The calculation in first and second and 
order in e is given in subsections 12.51 and 12.61 The general discussion in higher 
orders is presented in !2.7l The paper is closed by a remark on reparametrization 
with respect to the expansion in e (|2.8p . Although the presented arguments hold 
to arbitrary order in e, no solution in closed form is presented nor even a proof 
of convergency for the e-expansion. 



2 Equilibrium and Expansion 

2.1 Equilibrium 

Due to Archimedes' principle the body obeys 

V a = V h = ±V t , (6) 

where V a is the volume above, Vb the volume below the water plane and Vt the 
total volume. Condition $Q determines the height of the water plane. 
The potential energy V of the body is given by 

= pg(Z a -Z b ) = 2pg(Z a -^Z), (7) 

with m and p mass and density of the body, g the gravitational acceleration, 
z a and Zb the z-coordinates of the part of the body above and below the water 
plane, (the z-coordinate measured in direction of n, Z a , i?b, and Z t the integral 
of z over the volume above, below the water plane and over the total volume. 
In order that the body can float freely in all orientations the potential energy 
has to be independent of the orientation. Thus Z a — ^Z t has to be constant. 

2.2 Expansion 

In the following a unit vector u in d dimensions is decomposed into the com- 
ponent upward and a vector u' in the d — 1-dimensional space parallel to the 
water plane 

u = cos 9 n + sin 9 u' (8) 
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Then the element dild of the solid angle in d dimensions can be written 

dn u = dfVd6>(sm6l) d ~ 2 . (9) 

The angle G, at which the water plane intersects the surface of the body as a 
function of u' has to be determined from 

h = rcosQ. (10) 

Expanding 

Q = 1+SQ, r = r + Sr, 5r(Q) = ^&Lsr^ (11) 

(the Sr^ are the coefficients of a Taylor expansion of Sr in 50) one obtains the 
equation 

h = -(r + Y^^Sr^)sm5@, (12) 

3=0 J ■ 

which has to be solved for 50 for arbitrary u', which is not indicated explicitly. 
This eq. is invariant against the transformation h — > —h, 50 — > —50, — > 
(— y5r^\ thus (50 obeys the symmetry relation 

(50+ = -56_, (13) 
50+ := 5Q(h,{6r^}), 50_ := 56(-/i, {(-) J '5r^}). (14) 

The volume 14, and Z a are obtained by realizing that the volume, which is a 
hypersegment, can be divided into a hypersector (first integral) and a hypercone 
(second integral) 

K = i J dQ u ' J ( ] d0r d (9,u'){sin0) c 



\d-2 



h 



Ydr» u ,(r(0(u'),u')sin(0(u'))) d \ (15) 



(d- l)d. 

Z a = / dO u . / ( ) d0r d+1 (0,u')cos0(sin6') d - 2 

d+l J Jo 

- (d-l)(d + l) / d "»^ < 16 > 

It is practical to split the 0- integration into one from to 7r/2 and a second one 
from 7r/2 to 7r/2 + <50(u'). Then the expressions for V a and Z a are given by 

V, = VV+VM, Z a = Z^+Z^, (17) 
V™ = dn u , J*' 2 der d {e,u')(sm6) d - 2 , (18) 

Z« = -1- / dft u , / 7 d#r d+1 (0,u')cosfl(sin0) d - 2 ,(19) 
d + 1 J Jo 

= J dn u , g, Z^ = J dn u , g z (20) 
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with 



g(h,{Sr^}) = J ff -^-^( r (| + ( je + (u'))cos^e + (uO) d \ (21) 

I g (h,{6r if) }) = - AO' r d (- + 9' 1 {6r^})cos d ~ 2 9\ (22) 

d Jo 2 

g'(h,{8rW} = V- - (rf _ 1 ^ +1) (r(|+^e + (uO)cos( ( 59 + (uO) d ?2 ) 3) 



-1 



I r (/i,{<fr^}) = — L - / de'r rf+1 (^+^,{5rW})cos rf - 2 ^sin^. (24) 
a + 1 Jq 2 



Replace 



r (ZE + 0', {Sr ^}) = r (l - 6', {{-ySr^}) (25) 



and in a second step 9' by —6', Then I g reads 

I g = -=: / d0'r d (- + , ,{(-)^r«})cos d - 2 0'. (26) 

Insertion of (fT5|) yields for the integral 

I g (h,{6r^}) = -I g (-h,{(-ySr^}). (27) 

Similarly one shows 

I g *(h,{5rW})=I g ,(-h,{(-y8rW}). (28) 

Since 

r (~+6e + ,{8rW}) = r(~-6B+,{(-)t8rW}) = r(|+*e_, {(-)''<&•«>}) (29) 

one observes that the second terms in ([2~T|) and (j2"3"|) obey the same symmetry 
relation so that the relations 

g(h,{5r^}) = -g(-h,{(-) k 5r^}), (30) 
g*(h,{5r^}) = g*(-h,{(-y6r^}) (31) 

hold. In performing these calculations h depends on n, whereas 5r depends 
on u'. Since 5r is expanded around 8 — n/2 the Sr^ are to be taken at u' 
perpendicular to n. Thus g and g z depend on 

h = h(n), 5r (l) =Sr (J) {u'). (32) 

2.3 Expansion in ultraspherical harmonics 

The quantities under the f2 u ' integrals are expanded in ultraspherical harmonics 

oo 

/(u) = /(cos0n + S in0u') = £/;j(u), (33) 
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where the ultraspherical harmonics are eigenfunctions of the Laplacian on the 
unit sphere 

A u / ; ,(u) = -/(f + d-2)/ ; ,(u). (34) 

To distinguish this expansion from that in e the index / is always preceedcd by 
a semicolon. For given d and I there are 

#? = < 2, + d - 2 >T^2)f < 35 » 

linearly independent ultraspherical harmonics. 

The integral over O u / yields a function depending only on cos 9. The integral 
over / ; ; yields 

^{d/2-1), g\ 

dJW ;I (u) = EVi/ ;i (n p Md/2 T S / (36) 

(1) 

with the ultraspherical (Gegenbauer) polynomial 1 ^(x), which is the only 

function fa independent of u'. The factor in front is obtained, since the integral 
is the factor fld-i times the average over u', which for 9 — > approaches /;i(n). 
In appendix lA.3l a few formula for theses polynomials are listed. See also [TUl[TTj . 

(2) (2) 

In evaluating the integrals for Vw and Z\ (12"U)) one has to integrate over 
u' at 9 = f . Thus 

J dfi u ,/ ;I (u') = n d ^ d) f a (n), (37) 

l[ d) = % /2 1} (0) = H v^r(^) even 1 . (38) 
Cf /a-1) (l) 1 odd? 

These quantities are listed in appendix I A. 41 for d = 2 and 3. 

2.4 Vt, Z t , and Z« 

The total volume is given by 



dQ u < J d6r d (9,u')(im d - 2 9 



= >-xE(^(n)/ dg -' i^ ain^fl. (39) 
i 1/0 °; 

The integral over vanishes for all / except 1 = 0, 

V t = ^ d (r d ), (n). (40) 

There is only one ultraspherical harmonic Cq^ 2 *\x) — 1 for I = 0. {r d )-a is a 
scalar and does not depend on n. 
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Z t of the total volume is given by 

Z t = -r^—[dn u ,[ d9r d+1 (0,u')cos6sm d - 2 9 
a + 1 J Jo 



d8 cos e sin d - 2 9. (41) 



d+1 a -"Y 'J, C \ d/2 ~ 1] {\) 

The integral over 9 vanishes for all I except 1 = 1, 

* t = ^^) ;1 (n). (42) 

There are d ultraspherical harmonics for 1 = 1. They transform under rotations 
like the components of a vector. (r d+1 ) ; i(n) is the projection of volume times 
the vector to the centre of gravity onto n. 
One obtains 

r^(K W -^) = iE^(^(n) (43) 



with 

r(«) 



(d-2)Cf/ 2) (0) f even I 



and 



with 



^(^-^) = ^E^(^(n) (45) 

' (2Z + rf-2)C ; (d/2 - 1) (l) \l + d-l l-l J 

= I J 4v ^r(i±|±i) evenf . (46) 

[ odd I 

2.5 First Order in e 

Eq. (fl"!?|) yields in first order in e 

501 = (47) 



One obtains 

— = ^£// z V 1;i (n), (48) 
V,i=0, <£ = 0, Zg=0. (49) 
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Thus 

J; (2) ri ;i (n) =Ofor 1^0 (50) 

(z) 

has to be fulfilled. Since / ; v ' vanishes for odd I one can choose r\-\ arbitrarily 
for odd I, whereas r\-i = for even I with the exception of I = 0. ri ; o 7^ would 
change the radius. For simplicities' sake ri ; o = is chosen. 
Further one obtains 

1 1 T /^ _ J-iVrW, 



TO - > M ) = r^-^/Wn^n), (51) 



/„,, = -^-rj-', 9l = -^jr}- 1 , (52) 

^? = -3^-', (53) 

which yields 

h 1;l = (d-l)4 v) ri . tl . (54) 

Thus one obtains contributions hi-i proportional to r\-i for odd I, whereas h\-i 
vanishes for even I, since both r±- t i and vanish for even I. 

2.6 Second Order in e 

In second order in e one obtains 

SQ 2 = -^ + ^, (55) 



' {Z { 3-\z ta ) = rt^^^ + ^rl),, (56) 

1 



i 9 , 2 - -^rr^r, (57) 
a - -^-^rr, (58) 

(2) 1 ftfrg" 1 . (59) 



fi d -i a ' 2 2(d-l) 
This yields the equation 

jfW + ±rf), - ^^{hlh = 0. (60) 

For odd I one obtains 1^ — 0, (rf) ; / = 0, {h\).j = 0. Thus r 2 -i can be chosen 
arbitrarily, but in conformity with ([3]), the choice r 2 -i = is made for odd I. r 2 -i 
is determined for even I by 

T2-, = + (61) 
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One obtains for V 

L -(K ( ,2 - 5^,2) = 4-^l V) (r,r 2 + d —^rl), u (62) 



fi«j-i a ' z 2 



/*a = -^^^^(^-^ir^uO^iln)^- 2 , (63) 



.72 



2d(d 



1 T ,( 2 ) _ d+1 „ d _i, 



-K2 = ' iN rg-%(n) 



2d(d-l) 

- E^^^o" 2 ^)^^)^!^) (65) 
i 

F a 2 — §14,2 vanishes, since (ror2 + vanishes for odd Z and lj v ^ for even 

f 2) 

I. The second term in the expression for V a 2 vanishes, since (ri) ; ; vanishes for 

even Z and C; (d/2_1) (0) for odd I. Thus the condition (V^ + V$ - \V ta );i = 
reduces to /i2(n) = 0. 

Thus it has been shown that (fl]) and $2$ hold up to second order in e. Next 
it will be shown by complete induction that they hold in arbitrary order. 

2.7 Higher Orders in e 

gk and g^ are polynomials in h q and 5rp^ . They are linear combinations of 
terms 

nv.iWi m) ( 66 ) 

In order e fe contribute terms with 

fc = y^ffm +y^Pm- (67) 

Except for linear terms in the polynomials gk and {??, which have to be deter- 
mined, use of ([3J yields 

Ipr^>(-u') = (-)E„i»+E m ^ n^ m) ( u ')- (68) 

m m 

Thus only terms with even J2 m j™ + E m P™ contribute to and £, eq. 
(f2"0|) . From eqs. ([50)) and (f?T]) one deduces ^ q m + 3m 1S °dd for gk and even 
for gi. Due to ([3]) only h q with odd g can differ from zero. Thus all non-linear 
terms (n.l.t.) to V fc vanish for even k and such contributions to Z k ' vanish 
for odd k. 

Next consider the contributions — \V t ,k and — \Z t ^- Eq. (JTJ) yields 
(r d ) k (-u) = (-) k (r d )k(u), (r d+1 ) k (-u) = (-) k (r d +%(u). (69) 
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Since f-i(— u) = (— ) / ; i(u) there are only contributions (r d )k-i and (r d+1 )fc ; / with 
even k — I. 

Extracting the terms linear in and hk one obtains 

i 

I 

Consider odd k. Z^\ vanishes, the n.l.t. in Z^ vanish, too, since 1^ = 
for odd I and all n.l.t. contribute to odd I. One is free to choose r^i with odd 
I. One obtains hk-i from ([7T|l . Since all n.l.t. contribute to odd I, only hk-,i with 
odd I can differ from zero. 

Consider even k. Then all n.l.t. from Z contribute to even I . They determine 

(z) 

Tk-i with even I. For odd I the n.l.t. vanish and since D ' = for odd / one 
could choose non-zero r^i- However, choose r^i — for odd I and even k in 
agreement with ([3]). All n.l.t. to V ai k vanish. Moreover I^rk-i = 0, since the 
first factor vanishes for even I, the second for odd I. Thus hk = for even k. 
This completes the proof. 

2.8 Reparametrization 

In each order e k it is possible to choose r^i for odd I freely. We restricted 
ourselves to do this for odd k only. Basically it is sufficient to do this for k = 1 
only by reparametrizing 

oo 

^=E e2fer 2fc + l;Z ( 72 ) 

and r! 2k+1 . l = for k > 0, since it will yield the same result r(n). 

If one requires that the centre of gravity is located at the origin, then one 
has to adjust r^i in all odd orders k. Then in first order it has to be r 1;1 = 0. 
Thus the restriction ([3]) does not restrict the manifold of solutions. 

It is not necessary to introduce an expansion in e. Alternatively one expands 
r(n) + r(— n) — 2tq in powers of r(u) — r(— u). 



7=; [Z a --Z % )k - 

"d-1 ^ 



A Appendices 

A.l Auerbach's Solution 

Basically Auerbach[2] starts from the envelope of the water line. The length 21 
of the water line is constant. The line touches the envelope in the middle. Thus 
the envelope may be represented by 

x c (<f>) = .T + f d0's(0') cos(^'), VM = 2/o + / d^'s(^) sin(0')- (73) 
Jo Jo 
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The water line and thus the point, at which the envelope touches the water line 
has to be the same after <f> has increased by ir. Thus 

x c (ir) =x e {0), y e (w) = y e {0), s(w + <f>) = -sty) (74) 

is required, which yields 

x e (4> + n) = x e ((f)), y c ((f> + Tr) = y c {4>). (75) 

Then the boundary of the body is given by 

x(4>) = x c (4>) +£cos<j), y(4>) = y e (<j>)+£ sin <j>. (76) 

The scheme presented here is less elegant than the procedure used by Auerbach, 
however the present scheme applies to general dimensions. It would be desirable 
to have an elegant procedure which does not refer to a perturbation expansion 
for the construction of such bodies in d > 2 dimensions. 

A. 2 Envelope from h(n) 

For given h(n) one can determine the corresponding point r (n) on the envelope. 
It obeys 

(n + Sn) ■ r c (n) = h(n + Sn) (77) 
to first order in Sn, which is orthogonal to n. Expanding both sides 

n • r c + Sn ■ r c = h(n) + Sn ■ V/i(n) (78) 

yields 

r e (n) = nft(n) + Vft(n). (79) 

A. 3 Ultraspherical (Gegenbauer) Polynomials 

A few relations on ultraspherical (Gegenbauer) polynomials are given. Ro- 
drigues' formula 

r «( T ) - ( \n T(a + l/2)T{n + 2a) _ 2 , a+1/2 d" _ 2 y l+a -i/2 
°" W_l J 2-n\T(2a)T(n + a + l/2) {i X > dx« 1 ' 

(80) 

Recurrence relation 

2(n + a)xC:(x) = (n + l)C^ +1 (x) + (n + 2a- ^C^x). (81) 
Special values 

^(±i) = (±ir( n+2 ; _1 ). (82) 

C? m (0) = (-r( m+ 2~ i y C 2 Vi(0)=0. (83) 
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Integrals 



d-2 



yd-lMd/2) 



C^' (cos 6), l>0 



1(1 + d-2)" 

AO cos 0(sin 6) d - 2 c\ d,2 ~ 1] (cos 6)6.6 



The full solid angle is obtained from 



<>,, = <>,, i. / dflsin^ 2 ^^!-^ 1 ^ 2 ) 



r(|) 

with f2i = 2 for the two directions in one dimension, which yields 

n 2ix d l 2 

r(|) 

The duplication formula for the T-function reads 

r(2z) = ^2 2z - 1 T(z)T(z + h. 

\/7T 2 



(84) 



^(sin^- 1 , 1 = 6 



(86) 



(87) 



(88) 



A. 4 Dimensions two and three 
A.4.1 The limit d^2 

The normalization of C < f^ 2 ^ 1 \ ([50]) is not convenient for the limit d — > 2, since 
C approaches in this limit for / 7^ 0. One obtains the limit 

-,(d/2-l). 



lim 



cr'* - ^ (cos A) 



c, 



(d/2-1) 



(1) 



= T z (cos0) = cos(^), 



(89) 



where T; denotes the Chebyshev polynomials. That this normalization is not 
convenient, can also be seen from the factors (d— 2) in (|44l |46|) . However, one 
finds in this limit 



7, 



— 



r («) 
■/ 



even I 



odd I 



(-)' 


-)G-2)/2^i__ 



(-) 





(J-l)/2 1 
i 





(90) 
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A. 4. 2 Three Dimensions 



In three dimensions one obtains 



C; 1/2 (cos(9) = P z (cos 6>) 



(91) 



with the Legendre polynomials P\. One obtains 



even I 



odd I 





2'(/-l)(i+2)(i)! 2 



V. 




2i(i+i)!(i^i)! 







('-!)! 



(92) 
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